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Abstract. Thermodynamic properties and phase equilibria behavior for the hydrates 
and coexisting pairs of hydrates of common acids which exist in the stratosphere are 
assembled from new laboratory measurements and standard literature data. The 
analysis focuses upon solid-vapor and solid-solid-vapor equilibria at temperatures 
around 200 K and includes new calorimetric and vapor pressure data. Calculated 
partial pressures versus 1 IT slopes for the hydrates and coexisting hydrates agree well 
with experimental data where available. 


Introduction 

The depletion of ozone in the polar stratosphere is gener- 
ally accepted to be initiated by heterogeneous processes 
occurring on particulate surfaces [e.g., Solomon , 1988; Mo- 
lina , 1991, 1994]. These processes repartition chlorine from 
relatively inert and photostable “reservoir” compounds 
(principally HC1 and C10N0 2 ) to more photolabile forms 
(Cl 2 , HOC1) which are easily converted into the chlorine 
radicals (Cl, CIO), resulting in catalytic ozone destruction. 
Also key to extensive ozone destruction is the concomitant 
“denoxification,” that is, the repartitioning of nitrate species 
(NO x ) from the gas phase to the solid phase, where they are 
unavailable to decrease the chlorine radical concentrations 
by reforming the C10N0 2 reservoir species. 

Polar stratospheric clouds (PSCs) are broadly classified as 
being either type I, generally assumed to be crystals of nitric 
acid trihydrate (NAT), or type II, which are composed of 
ordinary ice [ Toon et al., 1986; Crutzen and Arnold , 1986; 
McElroy et al ., 1986]. On the other hand, the background 
“sulfate aerosol” of the midlatitude lower stratosphere 
consists of 60-80 wt % H 2 S0 4 aqueous solutions [e.g., 
Turco et aL , 1982]. Recent investigations have revealed that 
upon cooling to temperatures characteristic of the lower 
polar stratosphere (190-200 K), sulfuric acid droplets will 
swell not only from absorbing additional water vapor but 
also take up appreciable amounts of nitric acid from the gas 
phase and, upon continued cooling, become highly super- 
cooled sulfuric acid/nitric acid/water droplets [Zhang et al . , 
1993a] which crystallize at temperatures close to the ice frost 
point, eventually producing NAT and H 2 S0 4 hydrates [Mo- 
lina et al ., 1993; Beyer et al ., 1994], 

In addition to ice and NAT, liquid and frozen sulfuric acid 
solutions have been shown to efficiently catalyze chlorine 
activation reactions at sufficiently low stratospheric temper- 
atures [Molina et al 1993; Hanson and Ravishankara, 1993; 
Zhang et al., 1994]. Furthermore, liquid aerosol droplets 
participate in ozone destruction by enhancing the reparti- 
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tioning of oxidized nitrogen species via the heterogeneous 
reaction N 2 0 5 + H 2 0 —> HN0 3 , which occurs throughout 
most of the lower stratosphere [Kolb et al. , 1995]. 

Our knowledge of the physical properties and chemistry of 
these aqueous acids and their crystalline hydrates at low 
temperatures has greatly improved since the discovery of the 
Antarctic ozone hole. Vapor pressure measurements under 
conditions representative of the lower stratosphere have 
been carried out for the HN0 3 /H 2 0 [Hanson and Mauers- 
berger , 1988a, b; Worsnop et al., 1993], the HC1/H 2 0 
[Hanson and Mauersberger , 1990; Abbatt et al., 1992], and 
the H 2 S0 4 /H 2 0 binary systems [Zhang et al., 1993b]. Smith 
[1990] commented on the general thermodynamic consis- 
tency between HN0 3 /H 2 0 solid-vapor phase diagrams gen- 
erated from vapor pressure measurements and that expected 
from the standard thermodynamic properties of that binary 
system [Forsythe and Giauque , 1942]. In addition, thermo- 
dynamic data from various sources have been correlated and 
parameterized for the H 2 S0 4 /H 2 0 [e.g., Zeleznik , 1991] and 
HN0 3 /H 2 0 [Clegg and Brimblecombe , 1990] binary sys- 
tems. 

The acids considered here are H 2 S0 4 , HN0 3 , and HC1. 
On a molecule-per-molecule basis, bromine species are 
much more efficient in catalytic ozone destruction than those 
of chlorine, but atmospheric HBr mixing ratios are too small 
for significant amounts to be present in the condensed phase. 
The concentration of HF in the stratosphere is comparable 
to that of HC1; however, its hydrates are expected not to be 
stable under stratospheric conditions. They will likely have 
high HF vapor pressures as their stoichiometries are 
H 2 0 * nHF, with n > 1 [Cady and Hildebrand , 1930], rather 
than acid • aiH 2 0 as in the cases considered here. 

In this report we summarize phase equilibria of the indi- 
vidual hydrates as well as coexisting pairs of hydrates, with 
emphasis on the boundaries of the stability regimes, as 
predicted from literature thermodynamic data and in some 
instances augmented with new data from our laboratory. 
These include measurements of vapor pressures over coex- 
isting pairs of hydrates by mass spectrometry and determi- 
nations of the heats of fusion for the lesser studied hydrates 
HC1 • 6H 2 0 and HN0 3 * 2H 2 0 by differential scanning cal- 
orimetry. 
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Methodology for Solid- Vapor Equilibrium 

Consider the generalized binary process for two gaseous 
species forming a solid product C 


the hydrates. We begin with (1) for one hydrate and intro- 
duce 

v' a A -I- v’ b B <=> C' (6) 


v a A + v b B O C, (1) 

in which A and B are the gaseous components and the are 
integers that define the stoichiometry. The partial pressures 
of A and B , p A and p B , are related by the equilibrium 
constant K p for this process [e.g., Denbigh , 1981, section 
4.12] 


K P = ( 2 ) 

which is a function only of temperature. Second-order 
effects such as those arising from nonideal gas behavior or 
from the effect of the total pressure on the chemical potential 
of the solid are neglected in the derivation of (2). The 
equation is applied here only to the formation of binary 
hydrates (acid • «H 2 0); similar equations could be derived 
for three or more component hydrates. 

K p varies with temperature as 


for the other. As the equilibrium partial pressures of B over 
crystals C and C' are equal, the pair of equations, (2) for K p 
and the corresponding equation for K ' p , are combined to 
eliminate p B , which yields 



1 1 

T In Pa = — 1" Kp - ~T In K' , 
vb] v b v b 


with a temperature dependence 


(7) 


-RT*\ 

\ V B V B 


a i ; a \ d In p A 


dT 


= v a H° a - v b H%) 


/rzO t t/0 , uQx 

, (H c v A H a v bH b ) , 

V B 

or equivalently 


( 8 ) 


2 d In K p 

RT 2 

dT 


d In Kn A n 

= -R „ , x s + v b H?'$* s 

d \ T) A f ' A B J B 


- ^/;c ryStal = A// S ° ubl , (3) 

where R is the gas constant, the Hf s are the standard 
enthalpies of formation, and A//® ub | is the enthalpy change 
upon sublimation. The entropy changes, i.e., the intercepts 
of 


RlnK p = A 5 s ° ubl - AffVr (4) 

plots, may be found in a similar manner. 

Values for Hf s and 5/ s of the pure components at 298.15 
K and sometimes directly for the liquid solutions as well, 
such as for sulfuric acid, are given in standard compilations 
of thermodynamic data such as Lide [1993], Perry [1984], 
Chase et aL [1985], Dean [1992], etc. For other cases the 
heat of mixing is added to the sum of the individual heats for 
the pure liquids to obtain the total heat of formation of the 
mixture. As the temperature of interest here is not 298.15 K 
but near 200 K, adjustments are employed where heat 
capacity data, C p {T) y are available. Since only the difference 
between the heats of formation between the vapors and the 
solids is relevant, for convenience, only the components, 
mixtures, and crystals are adjusted to 200 K and not the 
elements from which they are formed; that is, rather than 
Hf s, we employ H f s at 200 K referenced to the elements still 
at 298.15 K. Thus, to obtain the heat of formation of a 
hydrate at 200 K, we consider the following cycle: mix the 
liquids at 298.15 K, cool to melting point, crystallize at the 
melting point, and further cool to reach 200 K 

V*200 K(cryst) = v a h};F* + v B Hfp M + 15 K H 

fTfusion r f 200 K 

+ Cp(T) dT + Ajjgyr + Cp(T) dT (5) 

J298.15 K JT tuiion 

In addition to the individual compounds (hydrates) we also 
consider the simultaneous equilibria of pairs of addition 
compounds (e.g., various hydrates of same molecule, differ- 
ing only in v f ), which further define the stability regimes of 


dl/T YA 

vb v ’b 

which represents the slope of a In p A versus 1 IT plot, that is, 
the equilibrium coexistence line for the two hydrates in 
question in a phase diagram. A similar equation can be 
derived for the intercepts (entropy changes). 

For the current example of water vapor over a pair of acid 
hydrates having n and n’ water molecules, we substitute 
H 2 0 = A, acid = B y v B = v' B = 1, v A = n , and v A — n' 
into (9): 


Atfsabl 


a//; 


subl 


V B 


V B 


(9) 


d In p h 2 o A//' subl - A H sM 
d(l/T) (n - n')R ' (1 ° 

Similarly, for acid vapor over the pair we have 

d In Padd n&H subl ~ n sub! 

dWT) (n' ~ ~n)R (U) 

Experimental Methods 

The general experimental methods and procedures used 
here have been described previously [Zhang et aL , 1993a, 
b], so only brief descriptions will be given here along with 
any details specific to this work. 

Vapor pressures were measured using a static technique in 
which the samples were held in a glass vessel with the 
temperature regulated in the range of 180-200 K. Partial 
pressures were determined with a quadrupole mass spec- 
trometer which had been calibrated for H 2 0 against pure ice 
and for HC1 against pressures extrapolated from the data of 
Fritz and Fuget [1956], For experiments involving the HC1 
hexahydrate, samples had to first be cooled to below 150 K, 
then warmed to induce nucleation. The estimated errors for 
these measurements are of the order of ±10%. 

Infrared transmission spectra of HN0 3 * 2H 2 0 were ob- 
tained by flattening a drop of the nitric acid solution (63.6 wt 
%) between silver chloride discs held in a temperature 
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Table 1 . Components at Standard State 


Species 

Hf 298.15 K, 
k J mol ' 1 , 
Liquid 

S 298.15 K. 
J mol' 1 K 
Liquid 

Hf 298.15 K, 
kJ mol “ 1 , 
Gas 

r200 K 
J 298.15 K 

Cp dT , 
kJ mol “ 1 , 
Gas 

H f 200 
K, u kJ 
mol “ 1 , 
Gas 

5 298.15 K, 
J mol -1 K 1 , 
Gas 

S 200 K, 

J mol -1 K _l , 
Gas 

h 2 o 

-285.830 ± .042 * 

69.950 ± .079'’ 

-241.826 ± .042'’ 

-3.395 r 

-245.221 

188.834 ± .042'’ 

175. 038 r 

HCl 

-167.15 ± A dr 

56.5 ± .2 dJ 

-92.31 ± .1 / 

-2.859‘ 

-95.169 

186.902 ± .005'’ 

175.271 r 

hno 3 

-174.18® 

155.679® 

-135.13® 

-4.687 r 

-139.817 

266.400® 

247.468 1 

h 2 so 4 

-813.989 ± .67* 

156.895 ± .08'’ 

-735.13 ± 8.4* 

-7.383 r 

-742.513 

298.796 ±2.1* 

269.003 r 


Error estimates are given where available for this and other tables. 
a Referenced to the elements at 298.15 K. 
b Chase et al. [1985]. 

c Using heat capacities from Chase et al. [1985]. 

J For dissolving HC1 gas in water to » dilution. 
e Lide [1993]. 
f Dean [1992]. 

9 Weast [1984]. 


controlled cell. Samples were prepared and loaded at room 
temperature, and spectra were taken at frequent intervals 
during the temperature cycle. A Nicolet 800 spectrometer 
with a MCT-A detector was used, scanning at 2 cm -1 
resolution. 

Enthalpies of fusion were measured by differential scan- 
ning calorimetry (DSC), using a Perkin-Elmer DSC-7. Tita- 
nium or gold plated capsules were used to hold the sample 
solutions. Calibration of the DSC was made using either ice 
or cyclohexane. We estimate the uncertainties in the en- 
thalpy measurements to be less than 10%. 

Sample solutions were prepared by diluting the concen- 
trated acids with deionized water. The compositions were 
analyzed by standard acid-base titration or density measure- 
ments at room temperature, with an estimated accuracy of 
±0.1 wt %. Typically, about 3 cm 3 of liquid were used in the 
vapor pressure measurements, while the sample volumes for 
the infrared and thermal analyses were less than 0.03 cm 3 . 

Results and Discussion 

Nitric acid hydrates. For this we use the Hf values given 
by hide [1993] and the heats of mixing the pure liquids at 
298.15 K given by Forsythe and Giauque [1942] as well as 
their heats of fusion and heat capacities. Note that care must 
be taken to consistently use values from the same source 
within a given calculation, as enthalpies of formation are 
much greater than those of vaporization, and hence small 
differences in Hf s would lead to large errors in A//s. The 
standard state values used here are summarized in Table 1. 

As mentioned above the temperature of interest here is 
around 200 K. Vapor phase enthalpies are adjusted from 
298.15 K to 200 K via the relation K H = Jj$ K C p 

dT using the heat capacity data given by Chase et al. [1985]. 
For the liquid solutions and crystalline hydrates we also use 
low-order polynomial fits to heat capacity data. On the basis 
of Forsythe and Giauque's [1942] heat capacity data, as 
illustrated in Figure 1, we calculate ^ p to be 

-11.367 kJ mol -1 forHN0 3 • H 2 Oii q and -13.909 kJ mol” 1 
for HN0 3 * 3H 2 0 liq . The value of C p dT is -2.997 kJ 
mol -1 for HN0 3 • H 2 O cryst and -8.464 for HN0 3 • 3H 2 O CITSt . 
Summing the contributions to (5), we obtain /// OOK = 
-505.45 kJ mol -1 for NAM and -1107.6 kJ mol -1 for NAT 
(referenced to the elements at 298.15 K). Table 2 was 


produced in this manner and Table 3 by a similar process for 
the entropies. 

For HN0 3 • 2H 2 0 we could not locate temperature de- 
pendent heat capacity data below 270 K, so there is greater 
uncertainty than for nitric acid monohydrate (NAM) and 
NAT. Liquid phase Cp values were extrapolated from 
Timmermans' [1960] compilation of binary systems proper- 
ties, and as a rough estimate for crystalline nitric acid 
dihydrate (NAD) we used a mean of NAM and NAT heat 
capacities over this temperature range. Using differential 
scanning calorimetry, we measured the enthalpy change in 
converting solid NAD at 232 K to liquid solution at 245 K to 
be -22.546 kJ mol” 1 from the baseline-corrected integral of 
the curve in Figure 2 from 232 to 245 K, that is, the combined 
heat of the peritectic melting of NAD to solution + NAT (the 



T(K) 


Figure 1. Heat capacity data used for liquid and crystalline 
HNOj * H^O (NAM), HN0 3 • 2H 2 0 (NAD), and 
HN0 3 * 3H 2 0 (NAT). Symbols for NAM and NAT mark 
data from Forsythe and Giauque [1942]. For the NAD melt, 
e.g., liquid ~65 wt % HN0 3 aqueous solution, data points 
are from Timmermans [1960], and, as an approximation for 
crystalline NAD, an intermediate value of NAM and NAT 
was used. 
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Table 2. Enthalpies 


mixing 

298.15 K From 


Species 

Standard 

Liquids 

Hf 298.15 K, 
Liquid 

^"fusion ’ 

°K 

iSea Cp dT . 

Liquid 

fusion at 
^fusion 

JflL ^ dT, 

Crystal 

Hj 200 K, 
Crystal* 

A TT 200 K 
subl 

h 2 o 


-285.83 ± .042'’ 

273.15 

-1.890" 

-6.011*' 

—2.420" 

-296.151 

-50.93 

HC1 • 3H 2 0 

+ 18.67 / 

-1006.0 

248.4 

-10.2 



- 1038.0 9 

-207.2" 

HC1 • 6H 2 0 

+ 9.48 / 

-1872.7 

203.2* 

-38.V 

—27.0' 

~0 

-1937.9 

-371.4 

hno 3 • h 2 o 

-13.56' 

-473.57 

235.5-' 

-ll.^ 

-17.5I-' 

-S-OtV 

-505.45 

-120.4 

HN0 3 • 2H 2 0 

—20.49^ 

-766.3 

238* 

-14.7' 

-20.15 m 

-4.4" 

-805.55 

-175.3 

HN0 3 • 3H z O 

-24.48-' 

-1056.2 

254.6' 

-13.9' 

-29. IE 

-8.46' 

-1107.6 

-232.2 

h 2 so 4 * h 2 o 

-27.85 c 

-1127.7 

281. 6 r 

-3.51" 

-19.93 ± .48 c 

-8.85* 

-1160.0 

-172.2 

H 2 S0 4 * 2H 2 0 

-43.72 c 

-1429.4 

233. 7 C 

-16.4° 

-18.08 ± .16 c 

-4.14 c 

-1468.0 

-235.0 

H 2 S0 4 • 3H 2 0 

-51.88° 

-1723.4 

236.8 f 

— 19.3* 

-23.37 ± .63 c 

-5.79* 

-1771.9 

-293.7 

H 2 S0 4 • 4H 2 0 

-57. \6 C 

-2014.5 

244. 9 C 

-20.3° 

-29.73 ± .9 C 

-8.17 c 

-2072.7 

-349.3 

H 2 S0 4 • 6»/2H 2 0 

-65.28 £ 

-2737.2 

220.3" 

-43.1" 

-33.07 ± 1.1 c 

-5.10° 

-2818.4 

-482.0 


Values are in kilojoules per mole. 

Referenced to the elements at 298.15 K. 
b Chase et at. [1985]. 

1 Zeleznik [1991]. 
d Lide [1993]. 

e Giauque and Stout [1936]. 
f Weast [1984]. 

9 Using Hanson and Mauersberger's [1990] measurements. 
h Vuillard [1957]. 

'This paper. 

i Forsythe and Giauque [1942]. 
k Metastable, extrapolated from Ji and Petit [1991]. 

1 Timmermans [I960]. 

m This paper, integral from 232 to 245 K minus the heat capacity integrals of NAD from 232 to 238 K and the melt from 238 to 245 K. 
n Using an intermediate value of NAT and NAM. 

0 Giauque et aL [I960]. 


peak beginning at 232 K) along with the heat of warming and 
melting the NAT + solution until no NAT remains (the peak 
ending at 245 K). 

Support for the presence of NAD in the calorimetry cell, 
rather than a NAT 4 - NAM mixture, comes from parallel 
experiments wherein samples of liquid HN0 3 • 2H 2 0 solu- 
tion were held between AgCl windows and cooled and 
rewarmed as described in the experimental section. As 
shown in Figure 3, the infrared spectra indicate the sponta- 
neous formation of NAD from the liquid upon cooling and 
the conversion of NAD to liquid + NAT on warming to the 
220-230 K range, followed by the complete melting of NAT 
by —245 K. Our spectra compared well with those reported 


for NAD by Ritzhaupt and Devlin [1991] and for NAD and 
“/3” - NAT by Koehler et aL [1992]. In addition, the bands 
in the 2400-2700 cm “ 1 region arising from the incorporation 
of -2% heavy water in the sample solutions aided in the 
identification of the NAD to NAT transition. In some cases, 
small amounts of NAT formed along with NAD upon cool- 
ing, and in these instances, of which Figure 3 is an example, 
the conversion of NAD to NAT began at lower temperatures 
than when no NAT formation was observed upon cooling. 

Though the actual process in the calorimeter was likely the 
conversion to NAT followed by its melting, as was observed 
in the infrared experiments, conversion of the measured 
^ 24 5 k H transition value to an approximation for the heat of 


Table 3. Entropies 


Species 

5° 298.15 K, 
Liquid 

f ^fusion 

J 298.15 K 

CplT dT, 
Liquid 

"n i 
if* 1 ' 4 

lZ*CplTd T. 
Crystal 

5 200 
K, 

Crystal 

Intercept of R In*, 
K p Versus MT at 
200 K, i.e., AS 

A SIR 

h 2 o 

69.950 ± .079* 

-6.6202 

-22.006 

-10.233 

31.091 

143.95 

17.313 

HCl • 3H,0 


-37.45 



136.0* 

564.4* 

67.88* 

HC1 • 6H 2 0 


-154.5 

-132.8 





hno 3 • h 2 o 

217. 00 c 

-42.78 

-74.38 

-13.75 

86.08 

336.4 

40.46 

HNOj • 2H 2 0 

280.2^ 

-55. 

-84.6 

-20. 

120.6 

477.0 

57.4 

HN0 3 • 3H 2 0 

347.15° 

-50.40 

-114.3 

-37.23 

145.2 

627.4 

75.46 

h 2 so 4 • h 2 o 

211.51* 

-12.15 

-70.77 

-36.82 

91.78 

352.3 

42.37 

H 2 S0 4 • 2H 2 0 

276.363* 

-61.85 

-77.37 

-19.09 

118.0 

501.3 

60.26 

H 2 S0 4 • 3H 2 0 

345.373* 

-72.58 

-98.71 

-26.53 

147.6 

646.5 

77.76 

H 2 S0 4 • 4H 2 0 

414.529* 

-74.97 

-121.4 

-36.75 

181.4 

787.8 

94.75 

H 2 S0 4 • 616 H 2 0 

587.819* 

-167.1 

-150.1 

-24.28 

246.3 

1160 

139.5 


Values are in J mol 1 K 1 . A S/R is unitless. Heat capacity integrals use data from same sources as in Table 2 
a Chase et al. [1985]. 

b Using Hanson and MauersbergePs [1990] measurements. 
c Forsythe and Giauque [1942]. 

d From parameterization of Clegg and Brimblecombe [1990], 
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-90 -70 -50 -30 -10 °C 



180 200 220 240 260 280 

T (K) 

Figure 2. Differential scanning calorimeter (DSC) trace for 
the conversion HN0 3 * 2H 2 0 to liquid upon warming. The 
first peak corresponds to the conversion of NAD to NAT + 
liquid, and the warmer peak is the melting of the NAT. The 
inset shows the temperature-composition phase diagram 
based on the work of Ji and Petit [1991], with dashed lines 
showing extensions of stable phase boundaries into metasta- 
ble regions. 


fusion of NAD at its metastable melting point at 238 K [. Ji 
and Petit , 1991] may be done by subtraction of the heat 
capacity integral of metastable NAD from 232 to 238 K and 
of the supercooled liquid from 238 to 245 K. This leaves 


-20. 145 kJ mol -1 for the (actual or hypothetical) melting of 
metastable NAD at 238 K. Equation (5) can then be applied 
to find Ay// 200 k * -805.55 kJ mol -1 . 

The slope of a In K p versus 1 IT plot for NAT at 200 K is 
calculated to be, using (3), (-232.2 kJ mol -1 )//? = -27930, 
which agrees very well with previously reported values. 
Worsnop et al . [1993] measured -28520 ± 470, Smith [1990] 
fit -28100 ± 500 to Hanson and Mauersbergers [1988b] 
measurements (and also calculated —27000 using a similar 
analysis to this one but without consideration of the heat 
capacities), and Mozurkewich [1993] calculated —27881 us- 
ing an analysis which included heat capacity data. For the 
intercept, A SIR, we find a value of 75.46 (with pressure in 
atmospheres), whereas Worsnop et al. measured 78.7 ± 2.7, 
Smith obtained 76.52 ± 1, and Mozurkewich 75.32. 

For nitric acid monohydrate we arrive at In K p = 40.46 - 
14480/7, which may be compared with In K p = (44.3 ± 1) 

- (15330 ± 500)/7 from Smith's [1990] analysis and In K p 
= (44.1 ± 1.7) - (15264 ± 300)/7 from Worsnop et al.' s 
[1993] observations. And for the dihydrate we obtain In K p 

- 57.4 - 21163/7, which can be compared to Worsnop et 
al.'s In K p = (64.9 ± 2.2) - (22630 ± 360)/7. 

The slope of In p Hz o versus 1/7 at the NAT-NAM 
coexistence, that is, the right hand side of (10), is calculated 
to be (-55860 J mol -1 )//?. For convenience, in comparison 
to a logjoP h 2 o = ^ plot, we have B = -55860/(/? In* 

10) = -2918, which agrees well with Hanson and Mauers- 
berger' s [1988b] measured value of -2819. Using (11), we 
estimate B = -3370 for the slope of a log 10 p H no 3 = A + 
BIT plot, which again agrees well with Hanson and Mauers- 
berger's [1988b] value of -3561. Other coexistence slopes 
are given in Table 4. 

Hydrates of hydrochloric acid. The table entries are cal- 
culated in the same manner as for the nitric acid hydrates. 



Figure 3. Infrared spectra taken through cycle of cooling and warming HN0 3 • 2H 2 0. This sample 
spontaneously crystallized to form the crystalline dihydrate upon cooling from 298 K to 203 K (traces a 
and b), and gradually converted to a trihydrate + liquid mixture upon warming to 227 K (traces c through 

f). 
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Table 4. Coexistence Line Slopes at 200 K 


Coexistence Pair 

Calculated Slope 
of log 10 PH ,0 
Versus l/T 

Experimental Slope 
of iog 10 p h o Versus 
l/T 

Calculated Slope 
of log JO P ac id 
Versus l/T 

Experiment Slope of 
log 10 Pacid Versus 

l/T 

ice:HCl • 6H 2 0 

-2660" 


-3440 

-3340* 

ice:HN0 3 • 3H.O 

-2660“ 


-4145 

— 3968 c 

ice:H 2 S0 4 * 4H 2 0 

-2660“ 


-7600 


ice:H 2 S0 4 • 6 '/i H z O 

-2660" 


-7880 


HCI • 3H 2 0:HCI • 6H 2 0 

-2860 


-2250 


HN0 3 • H 2 0:HNOj • 2H 2 0 

-2870 


-3420 


HN0 3 • H 2 0:HN0 3 • 3HiO 

-2918 

—2819° 

-3370 

-3561 r 

HN0 3 • 2H 2 0:HN0 3 • 3H 2 0 

-2970 


-3215 


H 2 S0 4 • H 2 0:H 2 S0 4 • 2HiO 

-3280 


-5720 


H 2 S0 4 • 2H 2 0:H 2 S0 4 • 4H 2 0 

-2985 

— 3236 rf 

-6310 


H 2 S0 4 • 4H 2 0:H 2 S0 4 • 6V2 H 2 0 

-2770 


-7155 



"Same as ice from assumption of pure crystal (Raoult’s Law). 
b This paper. 

c Hanson and Mauersberger [1988b]. 
d Zhang et al. [1993b]. 


The heat of vaporization of the trihydrate is estimated from 
the vapor pressure measurements of Hanson and Mauers- 
berger [1990]. The hexahydrate heat of fusion, -26.98 kJ 
mol -1 , was measured in our laboratory using differential 
scanning calorimetry, as for NAD. 

A plot of the HC1 partial pressure over the HC1/H 2 0 
system is presented in Figure 4. Vapor pressures over liquid 
solutions, extrapolated from Fritz and Fuget [1956], are 
shown as long dashed lines; the shorter dashed lines indicate 
where the vapor pressures of the liquids and solids are equal 
(ice, HC1 * 3H 2 0 or HC1 • 6H 2 0), i.e., the freezing enve- 
lopes; and the solid line defines the trihydrate-hexahydrate 
boundary based on measurements by Hanson and Mauers- 
berger [1990]. The triangles mark our measurements of p HCi 
over ice-liquid equilibrium mixtures, and the pluses are for 
ice-hexahydrate coexistence mixtures. We find excellent 
agreement between observed vapor pressures and those 
predicted using (11), shown as the dotted line in the figure, 
along the ice-hexahydrate coexistence line. 

Hydrates of sulfuric acid. Here the individual values for 
the K p are of less interest because solid-vapor equilibria are 
unlikely to be maintained at the extremely low ambient 
mixing ratios of H 2 S0 4 vapor, but the water vapor pressures 
over coexisting solids are of importance. As water pressures 
over ice-hydrate coexistence systems are negligibly different 
from those of ice alone (from the fact that the ice phase is 
very nearly pure), calculated and measured values only over 
the tetrahydrate-dihydrate coexistence mixture are pre- 
sented in Figure 5. 

Stratospheric Implications 

Phase diagrams presenting vapor pressures as a function 
of temperature (such as those shown in Figures 4 and 5) can 
be used to place constraints on the existence of various solid 
phases in the stratosphere. Consider a particular acid hy- 
drate: it is stable against evaporation only if its water and 
acid vapor pressures are equal or smaller than the corre- 
sponding environmental partial pressures. It was in this way 
that the formation of NAT clouds in the stratosphere at 
temperatures above the ice frost point was first predicted 
[Toon et al. , 1986; Crutzen and Arnold , 1986]. 

On the other hand, stability does not necessarily imply 


formation; nucleation needs to take place first if a condensed 
phase is to form from vapor or if a solid phase is to form from 
a liquid or from another solid. (Note, however, that under 
conditions such as those applicable to the atmosphere. 


Temperature (K) 



Figure 4. Plot of log P HCi versus 1000/T. As discussed in 
the text, the long dashed lines are for liquid solutions, with 
the concentrations in molal units labeled; the shorter dashed 
lines represent the freezing envelopes; and the solid line 
defines the trihydrate-hexahydrate boundary. Triangles 
mark our measurements of p HC] over ice-liquid equilibrium 
mixtures, and the pluses are for ice-hexahydrate coexistence 
mixtures. The dotted line is calculated from (11). 
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melting, evaporation, and sublimation always occur practi- 
cally at equilibrium.) This behavior leads to some nonequi- 
librium situations where phase diagrams are still useful, 
namely in terms of the stability of metastable phases. An 
example from Figure 4 is HC1 hexahydrate: it nucleates only 
from the liquid supercooled to below — 170 K [ Vuillard , 
1957; Abbatt et al., 1992]. Hence a practical phase diagram 
for a HC1/H 2 0 system that has not experienced tempera- 
tures below — 170 K (e.g., droplets in the stratosphere) is one 
in which the hexahydrate is absent; it contains regions in 
which the liquid, ice, or the trihydrate are metastable with 
respect to the hexahydrate, and yet valid predictions of 
melting, evaporation, and/or sublimation of the various 
phases can be made with a diagram that ignores the hexahy- 
drate, even if they are metastable with respect to the 
hexahydrate. Thus Figure 4 shows that unless some efficient 
nuclei for the hexahydrate are present, no HC1 hydrates are 
expected to form in the stratosphere, as HC1 partial pres- 
sures do not exceed 10 “ 6 torr and temperatures never fall to 
170 K. 

Metastability commonly arises in solid-solid phase trans- 
formations: if equilibrium is strictly maintained, as temper- 
ature and partial pressures change, a solid phase transforms 
into a second one when their coexistence line (i.e., as given 
by (10) and (11) for hydrates) is crossed. However, if nuclei 
for the second solid phase are absent and do not form 
readily, the first phase persists into the region of stability of 
the second one, until the first solid either evaporates or 
melts, depending on the partial pressure-temperature con- 
straints. Melting occurs under conditions represented by an 
extension of the freezing envelope of the first solid into the 
region of stability of the second one. 

We have been able to observe in the laboratory various 
types of metastable phase behavior, as expected from the 
phase diagrams. For example, ice melts in the presence of 
HC1 vapor at temperatures well below the eutectic (as 
illustrated in Figure 4 by the extension of the ice melting 
curve, the short dashes, past the intersection with the HC1 
trihydrate curve); neither the HC1 hexahydrate nor the 
trihydrate forms readily from ice [Abbatt et a/., 1992]. 
Similarly, when warming sulfuric acid tetrahydrate with the 
water partial pressure maintained at 4 x 10 ~ 4 torr, the 
tetrahydrate is not observed to transform into crystalline 
dihydrate along the H 2 S0 4 * 4H 2 0:H 2 S0 4 • 2H 2 0 coexist- 
ence curve at —217 K but is observed to become liquid at a 
slightly warmer temperature where the extended 
H 2 S0 4 * 4H 2 0:liquid curve is intersected [Zhang et aL, 
1993b; Middlebrook et al ., 1993]. On the other hand, if the 
two solid phases are initially present, it is possible to 
determine experimentally their coexistence line, for exam- 
ple, by changing temperature and simultaneously measuring 
the vapor pressures, provided, of course, that the timescale 
of the measurements is long enough for equilibrium to be 
maintained. It is in this fashion that the various coexistence 
lines described above were experimentally examined. Simi- 
larly, metastable phase boundaries are experimentally acces- 
sible by carrying out measurements in the absence of the 
appropriate nuclei and in the direction of equilibrium trans- 
formations, e.g., melting. 

Concluding Remarks 

The thermodynamic properties relating to vapor pressures 
of the acid/water binary systems important in the chemistry 


TEMPERATURE (K) 



iooo/kk') 

Figure 5. As in Figure 4 but forp H ^ 0 over the H 2 S0 4 /H 2 0 
system. The coexistence line is calculated from (10) and data 
points are from Zhang et al. [1993b]. 


of the stratosphere are now known fairly well, especially for 
the nitric and sulfuric acid systems. In addition to determi- 
nations of whether a given hydrate will be stable under 
certain conditions, accurate values for the equilibrium con- 
stants ATp’s are also necessary for calculating the rates of 
heterogeneous reactions occurring on the particle surfaces, 
as rates of reactions such as C10N0 2 + HC1 — > Cl 2 + HN0 3 
have been found to be highly dependent upon the water 
vapor pressure of the hydrate [Abbatt and Molina , 1992; 
Chu et al. y 1993; Hanson and Ravishankara , 1993]. 

We conclude from this analysis that the Kp's of the 
individual hydrates as well as the vapor pressures over 
coexisting pairs of hydrates which have been measured are 
in good agreement with values calculated using thermody- 
namic relationships. 
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